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Abstract 

A major drawback of most methods to find analytic expressions for soli- 
tary waves is the a priori restriction to a given class of expressions. To over- 
come this difficulty, we present a new method, applicable to a wide class of 
autonomous equations, which builds as an intermediate information the first or- 
der autonomous ODE satisfied by the solitary wave. We discuss its application 
to the cubic complex one-dimensional Ginzburg-Landau equation, and conclude 
to the elliptic nature of the yet unknown most general solitary wave. 

Keywords: solitary waves, cubic complex one-dimensional Ginzburg-Landau equation, 
Briot and Bouquet equations, elliptic function, genus, truncation. 
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1 Introduction 



When they are autonomous, nonlinear partial differential equations (PDEs) always ad- 
mit a reduction, called traveling wave reduction, to a nonlinear autonomous ordinary 
differential equation (ODE) defined in the simplest case by u(x, t) — U (£), £ = x — ct, 
with c a constant velocity. 

If the PDE is integrable, for instance by the method of the inverse spectral trans- 
form (1ST) £Q, then there exists a whole machinery to integrate the resulting ODE 
and find all the solitary waves (maybe with some requirement of decay at infinity). 
We will not consider this case. On the contrary, when the PDE is nonintegrable, the 
reduced ODE is also generically nonintegrable and serious difficulties arise, which we 
now describe on a much studied example. 

The one-dimensional cubic complex Ginzburg-Landau PDE (CGL3), 1 

iA t +pA xx + q\A\ 2 A-ijA = 0, pqj^O, (A,p,q) eC, 7 en, (1) 

which is usually written by physicists as an evolution equation in the case (Im p) (Im q) ^ 
0, 

t Q A t = rA + $(l + ia)A xx -(l+i(3)\A\ 2 A, (t ,r,£ ,a, 0) € K, (2) 

is a generic equation which describes many physical phenomena, such as the prop- 
agation of a signal in an optical fiber [2], spatiotemporal intermittency in spatially 
extended dissipative systems Ej ■ Its traveling wave reduction, 

A(x,t) ^ ^M{^)e l( -- ujt + i P(^\ i = x-ct, (c,cj,M,(p)eTZ, (3) 
defines the third order system in (M, ip'), 

M " M ' 2 /2 ( cM' \ . . . , , r n 
77 - ip -Si h 7 + s r (cip + uj) + d r M = 0, 

M' ( cM' \ 

ip" + iff— - S r ( — + 7 J - Si {op' + UJ) + diM = 0, 

with the notation for the six real parameters d r ,di,s r ,Si,g r ,gi, 

A _l -a 1 ■ 1 , • 1 + 1 22 

d r + idi = -, S r — ISt = -, Q r + IQi = h ~C S r . (5) 

p p P 4 

The mathematical solitary wave is the general solution, if it exists, of the system 
(QJ, equivalent to a single third order equation for M (see © below). 

In the CGL3 case properly said lm(p/q) ^ 0, the system (0J is only partially inte- 
grable, which means that the general solution (which should depend on three arbitrary 
integration constants) does not exist, and the question is to find, in closed form, the 
solution which depends on the largest possible number of arbitrary constants, let us 
call it the general analytic solution. This number of arbitrary constants, smaller than 
three, can be computed either from singularity analysis @J[i], or from topological 
arguments and it is equal to one, namely the origin £o of £. The question is 
therefore to find the unique particular solution of pjl. 

M = /(£-&), (6) 

which has codimension zero (the codimension is defined as the number of constraints 
among the parameters dr,di,s ri Si,^,c,uj of the system @). 

At present time, all existing methods (mainly the so called truncation methods, 
see e.g. the summer school lecture notes ^S|) have failed to find this codimension-zero 
solution, they have only found several solutions with a nonzero codimension in which 

1 Only the imaginary part of 7 can be absorbed in the definition of A, not the real part. As to the 
contribution of a group velocity term —ivA x , it can be absorbed in the definition of t since v is real. 
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M is a polynomial in a variable which satisfies a Riccati equation. What we present 
here is not just an additional method which could maybe yield a new particular 
solution with a nonzero codimension, but an entirely new approach which for sure 
will find, sooner or later, the desired expression /(£ — &>) m ©• Rather than looking 
for an explicit, closed form expression, we look for the first order ODE resulting from 
the elimination of £o between (|HJ| and its derivative. Because of the invariance by 
translation of £, eliminating £o is equivalent to eliminate £ — £o, therefore this first 
order ODE is autonomous. There exist classical results of complex analysis, mainly 
due to Briot and Bouquet ^3 pages 58-59], which establish the following important 
points: 

1. the knowledge of the above mentioned nonlinear first order autonomous ODE 
implies the knowledge of the expression /(£ — £o)j 

2. the nature of this solution /(£— £o) can only be elliptic (i.e. meromorphic doubly 
periodic in the complex plane) or a degeneracy of elliptic (i.e. trigonometric 
(simply periodic) or rational). 

Therefore, if one is able to directly build this first order ODE, this suppresses the 
need for making a good guess for the class of expressions /(£ — £o), since this class is 
now an output, not an input, of the method. 

The new method presented here is particularly suited to nonintegrable equations, 
it is algorithmic, and is very easy to implement since it only requires the resolution 
of a linear system. 

The paper is organized as follows. 

In sectional we recall the already known particular solutions of (JHJ) in the CGL3 
case, together with their nonzero codimension. In section|3l we give the corresponding 
first order equations which they satisfy. In section 0] we briefly recall the previous 
methods, then we summarize the relevant classical results. In sectional we present the 
new method and its advantages. In section® we give the results for CGL3. Finally, 
in section we indicate the class of dynamical systems to which the method can be 
applied. 



2 Particular solutions with a nonzero codimension 

Rather than the system of two coupled ODEs (@J , it is necessary, in order to apply our 
method, to consider a single ODE, e.g. for the field M = \A\ 2 , obtained by eliminating 
<p between the system (@J, 



, cs r G' — 2cSjG ( , cs r \ 2 G 
9 = 



'_ ( I _ Cfr\ 2 = G_ 

2 ' 2M 2 (g r -diMY V 2 J M 2 ' 1 ' 

(G 1 - 2cs t G) 2 - AGM 2 {d,M - g r ) 2 = 0, (8) 

1 1 cs- 

G = -MM" - -M' 2 - -^MM' + d r M 3 + g t M 2 . (9) 

In the CGL3 case properly said lm(p/q) ^ 0, only three solutions are currently 
known. In all of them M is a second degree polynomial in r = (fc/2) tanh/c£/2, with 
k 2 a real constant. These are 

1. a codimension-one source or propagating hole P] (heteroclinic orbit), 

2. a codimension-two pulse or solitary wave |17| (homoclinic orbit), 

3. a codimension-two front or shock |15j (heteroclinic orbit). 

To avoid carrying heavy expressions, let us make the following nonrestrictive sim- 
plification. Out of the five parameters d r ,di,g r ,gi,csi of the ODE JHJ), only three 
are essential (g r ,gi,c, equivalent to 7, cj,c). Indeed, just like for NLS, p and q (i.e. 
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d r + idi and s r — is,) can be rescaled to convenient numerical values, making rational 
all the coefficients of the above solutions, such as 



p = — 1 — 3i, q = 4 — 3i, 

1 , 3 1 3 

d r = ~, di = -, Sr = -J 6 , Si = --. (10) 

Each of the three solutions depends on an additional sign and, for these numerical 
values, the two sets of three solutions are 
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Remark. It has been predicted [S] the existence of a fourth physically interest- 
ing solution, which is a codimension-one homoclinic hole solution with an arbitrary 
velocity c. From the study of the structure of singularities of the CGL3 equation 
[S] , the most likely analytic expression to represent it would be a polynomial in both 
tanhfc£/2 and sechfc£/2 of global degree two, 



M 



3 V&ff + 8d 2 
2d? 



9d r 

~2dl 



tanh+ci | sech+^-^ tanh 2 +C3 tanh+C4, (17) 



with the homoclinic condition C3 = 0. We have checked that, unfortunately, such a 
solution does not exist in the CGL3 case Im(p/g) ^ 0. 



3 The corresponding first order subequations 

The only arbitrary constant in the two sets of three solutions of CGL3 is the origin £0 
of £, which represents the location of the movable pole (movable means which depends 
on the arbitrary initial conditions). By eliminating it between the expressions for M 
and M', one obtains first order ODEs. For the set — (fl3|) . these are 

( M ' + l M+ +2 ( M + b) ( m - si - K> 2 = °> (18) 



4 



M + 2(M - g r )M 2 = 0, (19) 
(M' + |m) 2 + 2M 3 = 0, (20) 

with the respective constraints on (g r >9i>c) already indicated, and similarly for the 
set jg-JTSJ, 

M '+ C -M-^] - (m--^} (M+-^-*g r ) = 0, (21) 

/2 



5 500 / V 100 / V 100 



2 
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M'-Mi^M--g r j =0, (22) 

M' + -mY - M 3 = 0. (23) 
5 / 

These six first order ODEs will be called subequations of the third order ODE (JSJ, 
since (JHJ is a differential consequence of any of the six first order ODEs, with the 
respective constraints on (<7 r , g^, c). 

Another information, quite useful to get some hint on the sought expression /(£ — 
Co) m ®, is made of the similar expressions for the unique integrable case of the 
CGL3 PDE QJ, namely q/p e R,7 = 0, which defines the nonlinear Schrodinger 
equation (NLS), 

%A t + pA xx + q\A\ 2 A = 0, pq ^ 0, A e C, (p, q) e K. (24) 

In this case, as opposed to the generic CGL3 case, the general solution of the system 
(0} exists. Indeed, this system admits two first integrals K±,K2, 

Mip'-yM = K 1 . (25) 

so the general solution is the singlevalued closed form expression 

M=-2^(p(£-£o,S 2 ,S3)-eo), <ff - £ = (27) 
q 2p M 

eo = K\ = - (E) 2 (4e3 - 92 e - 93 ), K 2 = -d ( 6 e* - f ) , (28) 

in which p is the Weierstrass elliptic function defined by 

p' 2 = Ap 3 - g 2 p - g 3 . (29) 

It depends on four fixed constants, p, q, c, u> (fixed means which appears in the defi- 
nition of the ODE), and three movable constants. Out of these three, one (£o) just 
represents the translational invariance, and the two important ones are (32,93) or 
equivalently (Ki,K 2 ). 

The two physically meaningful solitary waves of NLS (the bright one and the dark 
one) are obtained from this mathematical solitary wave by imposing some decaying 
conditions at infinity on the real axis of the complex variable £. The same feature 
will certainly apply after the sought solution of CGL3 has been found. 

4 The classical results on first order equations 

The six subequations of CGL3 all have the same degree two (in the highest derivative 
M'), just like the (unique) similar subequation (|26|l for NLS, but the difference with 
NLS is the absence of first integrals K\, K 2 . 
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The six subequations are by construction in a one-to-one correspondence with 
the six solutions of section [21 therefore it is equivalent, in principle, to search for 
subequations or for solutions. However, there is a big difference between these two 
classes of methods. 

In every method which directly searches for solutions, one must a priori assume a 
given class of possible solutions, then compute whether there indeed is some solution 
in this class. All these methods can therefore be qualified as "sufficient" . The main 
existing such methods are 

1. the one-family truncation method [221 El which assumes the class of polynomials 
in tanh, here 

M = c 2 t 2 + cir + c , c 2 ^ 0, t = (k/2) tanh fcf/2, (30) 

2. the two-family truncation method f) ;i App. A] ^3], which assumes the class of 
polynomials in tanh and sech, here 

M = (d lT + d Q )cr + c 2 t 2 + Ci t + c , (d u c 2 ) ^ (0,0), (31) 
t = (fc/2)tanhfcf/2, a = {k/2) scchfc£/2, (32) 

3. the elliptic method ^j, which assumes the class of polynomials in the 
Weierstrass function p(£) and its derivative p'(£), here 

M = Cl p+c , ci ^0. (33) 

Therefore it is easy to build an ODE which escapes these methods, e.g. any dif- 
ferential consequence of the first order ODE 

[VI ' 2 i (io?i#2 3 ) flr / , Qfi „ # 4 17 ji/r2 1 

whose solution is 



M + [ 12M 2 --\M' + 36M 4 - — M 2 + ^ = °> (34) 



M = tanh «-&) . (35) 
2 + tanh 2 (£-£o) 

On the contrary, the search for first order autonomous subequations requires no a 
priori assumption at all, and, from classical results which we now recall, the knowledge 
of the first order subequation is indeed equivalent to the knowledge of the explicit 
expression ©. As opposed to the previous methods, which are "sufficient" as said 
above, the proposed method can be qualified as "necessary" . 

Consider the class of first order nonlinear algebraic ODEs, 

F(M, M') = 0, (36) 

in which F is a polynomial of two variables. The unknown first order ODE resulting 
from the elimination of £o between © and its derivative, under the (reasonable) 
assumption that it is algebraic, belongs to this class. 

A first theorem was established by Briot and Bouquet (see, e.g., Vol. II, §139 
p. 284]): If the general solution of the autonomous ODE ]86}) is singlevalued, then 
this solution is either an elliptic function, or a rational function of e ax , a being some 
constant, or a rational function of x. Note that the third case is a degeneracy of the 
second one, itself a degeneracy of the first one. 

A second result, of immediate practical use, is due to Painleve pages 58-59]. // 
the general solution of the autonomous ODE \Sb)) is singlevalued, then the necessary 
form of this ODE is 

m 2m — 2k 

F(M, M') =J2 Yl = 0. a o,™ = 1, (37) 

fc=o j=o 

in which m is a positive integer and the aj^k's are constants. 

We now present a method which, given the integer m, directly computes the 
coefficients of <{3"?|) . 
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5 The linear method based on the Laurent series 



The general analytic solution of the ODE JBJ), which depends on the single movable 
constant £o , is known locally, this is the Laurent series E] 

M± = ^|^X- 2 (l + f X + 0( X 2 )) , X = £ - 6>, (38) 

in which £o is the location of a movable double pole, and A is the positive fixed 
constant 



A = ^/9d 2 + 8dl (39) 

If one could eliminate £o> i- e - £ — £cb between the Laurent series and its deriva- 
tive, this would settle the question of establishing the codimension-zero first order 
subequation. 

The only piece of information we have on the unknown global expression /(£— £o) is 
its local representation by this Laurent series. In particular, there is no mathematical 
proof yet that /(£ — £o) is globally singlevalued. However, a numerical investigation 
[25] of the singularities of / in the complex plane displays a pattern of singularities 
which looks like a doubly periodic set of poles and zeros, i.e. the signature of some 
elliptic function. Let us therefore assume the singlevaluedness of /. 

The algorithm to convert this local information (the Laurent series) into a global 
one (the codimension-zero first order subequation) is the following. 

1. Choose a positive integer m and define the Briot and Bouquet first order ODE 
l[37[l. It contains (m + l) 2 — 1 unknown constants aj,k- 

2. Compute J terms of the Laurent series, with J slightly greater than (m + 1) 2 — 1, 

(40) 




with in the CGL3 case p = —2. 

3. Require the Laurent series to satisfy the Briot and Bouquet ODE, i.e. require 
the identical vanishing of the Laurent series for the lhs F(M, M') up to the 
order J 



F ee X D V F jX j + 0( X J+1 ) , D = m(p- 1), Vj : F j = 0. (41) 




If it has no solution for o^fe, increase m and return to first step. 
4. For every solution, integrate the first order autonomous ODE (|37|l . 

The core of the method is the third step, in which the system of J + 1 equations 
Fj = in the (m+ 1) 2 — 1 unknowns is linear and overdetermined, therefore quite 
easy to solve. Indeed, it has generically no solution until m reaches its correct value. 

As to the fourth step, it is also quite easy. Indeed, when the first order autonomous 
ODE (|S7|l has the property that its general solution is singlevalued, there exists a clas- 
sical method due to Poincare to perform the integration, and it has been implemented 
at least in one computer algebra language One first has to compute the genus 
of the algebraic curve F(M,M') = 0, which can only be zero or one because of the 
singlevaluedness. If the genus is zero, the general solution of Q37JI is a rational function 
of e a ^, with a constant, or a rational function of £. If the genus is one, the general 
solution is an elliptic function of £, which can always be put in the canonical form of 
a rational function of p(£) and p'(£). 
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The above method clearly includes the three previous ones mentioned in Section 
13 which are recovered with the respective values m = —p, — 2jj, —2p. Therefore it can 
only find more results, and more easily since the main computation (third step) is 
just linear. 

Remark 1. In addition to the Briot and Bouquet selection rule j + 2k < 2m, one 
also has the selection rule m(p — 1) < jp + k(p — 1), which, when p ^ —1, further 
restricts the number of coefficients &j,fc. 

Remark 2. There is no upper bound on m, and this is the only inconvenient of 
the method. However, there exists a lower bound, obtained by the requirement that 
the most singular term M' m ~ M™x m<i>-1 ' ) must balance some other term in <|37[) . 
For p — — 2 (the CGL3 or NLS case), this forbids m = 1 (because M' cannot balance 
M 2 , nor M, nor 1), and the loop must start at m = 2. 

6 Results for CGL3 

With the numerical values (|10|l . one has A = 9/2, and the two Laurent series are 

"-=^(- 2 +Mf-f-^)* 2+o( * S) )' (42) 

The existence of two Laurent series, rather than just one, is a feature which the 
subequation must also possess, and this has the effect of setting the lower bound to 
m = 4 instead of 2. Indeed, the lowest degree subequations 

F 2 = M' 2 + M'(ai,iM + a ,i) + a 3 , Af 3 + a 2 , M 2 + ai fi M + a ,o = 0, (44) 
F 3 = M' 3 + M' 2 (a h2 M + o , 2 ) + M'(a 3 ,iM 3 + a 2>1 M 2 + a M M + a ,i) 

+a 4 , M 4 + a 3 , M 3 + a 2 , M 2 + a lfi M + a 0:0 = 0, (45) 

have the respective dominant terms M' 2 +a3$M 3 and M' 3 + a^^M 1 'M 3 , which define 
only one family of movable double poles. 

To give a detailed account of the method, let us nevertheless start with the lowest 
bound 771 — 2, for which l|44(l can only be satisfied by one series, e.g. I|42l) . thus 
preventing the full desired result to be obtained. The six coefficients aj k of (|44|l are 
first computed as the unique solution of the linear system of six equations Fj = 0, j = 
0,1,2,3,4,6. Then the J + 1 — 6 remaining equations Fj = 0,j = 5,7 : J, which 
only depend on the fixed parameters (g ri gi,c), have the greatest common divisor 
(gcd) 3gi + 2g r + 3c 2 /50, and this factor defines the first solution (fTTfl . namely the 
codimension-one propagating hole of Bekki and Nozaki. After division par this gcd, 
the system of three equations Fj — 0,j — 5,7,8, provides two and only two other 
solutions, which are flt^l and pHJl. with the respective constraints (c = 0,gi = 0) and 
(g r = 0,50<7i — c 2 = 0), and all the remaining equations Fj — 0,j> 9, are identically 
satisfied. 

Therefore, with this lower bound m — 2, one already recovers all the presently 
known first order subequations. Finally, for each of the three subequations, the fourth 
step finds a zero value for the genus and returns the general solution as a rational 
function of e a ^~t°\ which basic trigonometric identities then allow to convert to the 
second degree polynomials in (fc/2) tanhfc(£ — £o)/2 listed in ((TT|l - l(r^|) . 

With the correct two-family lower bound m = 4, which corresponds to 18 un- 
knowns dj t k and at least 24 terms in the series, we have checked that there is no 
solution other than the above three. This situation is quite similar to the absence of 
solution in the class lfTT|l . and it just reflects the difficulty of the CGL3 equation. 

The case m = 8 (60 unknowns a^fe and at least 66 terms in the series) is currently 
under investigation but preliminary results seem to indicate the absence of any new 
solution, and we are now automatizing the computer algebra program in order to 
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handle much larger values of m. Indeed, if one admits the singlevaluedness of /(£ — £o) 
in |JBj, then for sure there exists some positive integer m , alas without upper bound, 
at which the result will be obtained. 

7 Domain of applicability of the method 

As we have seen, the present method contains the three main methods already avail- 
able, and its cost is minimal since the main step is a linear computation. 
The two key assumptions behind this "subequation method" are, 

1. a Laurent series should exist, 

2. a first order autonomous algebraic subequation should exist. 

Its best applicability is therefore nonintegrable N-th order autonomous nonlinear 
ODEs admitting a Laurent series which only depends on one movable constant, such 
as the present CGL3 ODE (JSJ or the traveling wave reduction of the Kuramoto- 
Sivashinsky equation ESj ■ 

Two examples of inapplicability are 

1. the Lorenz model, in which the Laurent series generically does not exist and 
has to be replaced by a psi-series |20) . 

2. the autonomous ODE M'" - 12 MM' -1 = 0, which admits the first Painleve 
transcendent as its general solution, a case in which no first order subequation 
exists. 

8 Conclusion 

The advantage of the Laurent series is to exclude the contribution of chaos in a 
dynamical system, its disadvantage is to be a local piece of information. The present 
method converts this local information into the global one of a first order autonomous 
subequation, without any other assumption than the singlevaluedness of the general 
solution of this first order ODE. This is therefore a significant improvement over 
previous methods, which all require an additional assumption restricting the sought 
solutions to some class of analytic expressions. 

From classical results, and again under the single assumption of singlevaluedness, 
we conclude that the yet unknown codimension-zero solution is elliptic, i.e. equal to 
some rational function of p(£ — £o) an d — £o)- 
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